Low-momentum interactions in three- and four-nucleon scattering 
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Low momentum two-nucleon interactions obtained with the renormalization group method and 
the similarity renormalization group method are used to study the cutoff dependence of low energy 
3JV and 4JV scattering observables. The residual cutoff dependence arises from omitted short-ranged 
3iV (and higher) forces that are induced by the renormalization group transformations, and may help 
to estimate the sensitivity of various 3N and 4iV scattering observables to short-ranged many-body 
forces. 
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I. INTRODUCTION 

Modern few- and many-body calculations of nuclear 
structure and reactions are based on the picture of point- 
like nucleons interacting via two- and three-nucleon po- 
tentials. For this purpose, a number of high precision 
(% 2 /datum ~ 1) but phenomenological meson exchange 
models of the two-nucleon (2N) force such as the Ni- 
jmegen Q, Argonne V18 (AV18) and CD-Bonn Q 
potentials have been developed over the past decade. 
However, with phenomenological models it is not clear 
how to construct consistent three-nucleon (3A) forces 
and other operators. The lack of a systematic organi- 
zation or counting scheme results in model-dependent 
predictions, as there is no way to make controlled com- 
parisons between the different force models. More re- 
cently, substantial progress has been made in construct- 
ing nuclear interactions from chiral effective field theory 
(EFT) d, which is based on the most general local 
Lagrangian with nucleon and pion fields and all possible 
interactions consistent with the (broken) chiral symme- 
try of QCD. In contrast to phenomenological interaction 
models, the EFT approach is universal and provides a 
model-independent framework with a systematic organi- 
zation of consistent 2N, 3N and higher-body forces (and 
other operators) prescribed by the power counting. 

For both phenomenological and EFT potentials, nu- 
clear few- and many-body calculations are complicated 
by strong short-range repulsion and tensor forces that 
necessitate highly correlated trial wave functions, non- 
perturbative resummations, and slowly convergent ba- 
sis expansions. However, the non-perturbative nature 
of inter-nucleon interactions is strongly scale dependent 
and can be radically softened by using the renormaliza- 
tion group (RG) to lower the momentum cutoff that is 
present in all nuclear interactions. A consequence is that 
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many-body calculations become much more tractable 
at lower resolutions, resulting in calculations that are 
amenable to straightforward perturbative methods, sim- 
ple variational ansatze, and rapidly convergent basis ex- 
pansions [IjHBUl- The RG approach has the important 
advantage of being able to vary the cutoff as a tool to op- 
timize and probe the quality of the many-body solution, 
and to provide estimates of omitted terms in the Hamil- 
tonian. 

The above considerations have motivated the con- 
struction of low-momentum potentials Vlnw k through 
the renormalization group (RG) method [10, E| and, 
more recently, by the similarity renormalization group 
(SRG) method gQ. Both methods serve to eliminate 
the strong coupling between low- and high-momentum 
modes in the Hamiltonian such that low-energy observ- 
ables are preserved. In the RG method, one integrates 
out the problematic high-momentum components of the 
input interaction above a momentum cutoff A, leading to 
a new energy independent potential V\ OV! k that has the 
same low-energy on-shell transition matrix (t-matrix) as 
the input potential. In the original approach A consti- 
tutes a sharp cutoff above which the t-matrix is zero; the 
method has since been generalized to include a smooth 
momentum-space regulator to avoid technical difficulties 
stemming from the sharp cutoff [8J. On the other hand, 
the SRG method uses a continuous sequence of unitary 
transformations that weakens off-diagonal matrix ele- 
ments, driving the Hamiltonian towards a band-diagonal 
form [a,0|- In contrast to the RG method, SRG preserves 
both low- and high-energy observables independent of the 
value of the flow parameter A that provides a measure of 
the spread of off-diagonal strength. However, as with 
the standard RG, the calculation of low-energy observ- 
ables is decoupled from the high-momentum physics with 
SRG-evolved potentials (i.e., one can truncate intermedi- 
ate state summations to low momenta without distorting 
low-energy observables). 

Observables are scale-independent quantities. It is 
well-known that RG (SRG) transformations generate 
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short-range many-body forces (in principle, up to A- 
body) that "run" with the cutoff to maintain exact A (A) 
independence of A-body observables. If the RG transfor- 
mation is truncated at the 27V level, then the resulting 
cutoff-dependence in 37V observables may provide an es- 
timate of omitted short-range 37V forces in the Hamilto- 
nian. Along these lines, low-momentum 27V potentials 
have been recently used in three- and four-nucleon(47V) 
bound state calculations [l2j as a means to assess the 
size of omitted higher-body forces by varying the cutoff. 
There, it was found that the induced 37V forces due to the 
truncation to low momentum are of the same order as the 
so-called "bare" 37V forces attributed to integrating out 
excitations of nucleons. That is, the cutoff-dependence of 
the 37V binding energies was rather weak, varying by only 
1 MeV over a large cutoff range, which is comparable to 
the 0.7-1 MeV binding provided by the missing "bare" 
37V forces in conventional models and EFT calculations. 
In this sense, the RG evolution to low momentum does 
not induce strong short-ranged three-body force contri- 
butions to these 37V bound state observables. Similar 
results were obtained in 47V bound state calculations, 
where the various 27V Vi ow k calculations did not differ 
any more from the phenomenological Tjon-line than did 
calculations using 27V plus adjusted 37V forces. 

In the current study, we extend the cutoff-dependence 
study of Ref. [l2| to 37V and 47V scattering observables. 
In particular, we apply RG- and SRG-evolved 27V inter- 
actions to study how the neutron-deuteron (n-d) elastic 
vector analyzing power A y and the space star cross sec- 
tion in n-d breakup change with the cutoff. These being 
the two major long-standing failures of realistic interac- 
tions in their description of 37V data at low energy, one 
would like to use cutoff-dependence as a tool to assess the 
sensitivity of these observables to omitted short-range 37V 
force effects. Likewise, the same applies to observables in 
47V scattering that show large deviations to data, namely 
the total neutron-triton (n-t) cross section cr t around the 
resonance region at neutron lab energy E n — 3.5 MeV 
and the p- 3 He A y that also misses the data by as much 
as 25 - 40%. 

In Section |TT] we study 37V observables and in Sec- 
tion [TTT| 47V observables. Finally in Section ITVl we present 
the conclusions. 



II. THREE-NUCLEON OBSERVABLES 

The results shown in this section are obtained from 
the solution of the symmetrized Alt, Grassberger and 
Sandhas (AGS) equations [HI for the 3N system using 
the numerical techniques of Ref. [T3 |. In order to re- 
late the present work to the findings of Ref. [12| we re- 
peat in Fig. Q] the cutoff dependence of the triton bind- 
ing energy e t for CD-Bonn, AV18 and EFT potential at 
next-to-next-to-next-to leading order (N3LO) [f| based 
Mow k potentials using RG (left side) and SRG (right 
side) methodologies. In contrast to the calculations of 
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FIG. 1: (Color online) Triton binding energy as function of 
RG cutoff A (left side) and SRG parameter A (right side). 
Results derived from CD Bonn (solid curves), AV18 (dashed 
curves), and N3LO (dotted curves) potentials are shown. The 
horizontal line at et = —8.482 MeV is the experimental value. 



Ref. [l 21 ] with a sharp cutoff A, for simpler numerics we 
use a smooth regulator of the form exp (— (/c 2 /A 2 )f). The 
results are consistent with the ones of Refs. (5.1l2l|. At 
first glance, the SRG parameter A that provides a mea- 
sure of the spread of off-diagonal strength is not obvi- 
ously related to the cutoff A in the RG. However, in 
Ref. [7J] it was found that the "decoupling scale" for SRG- 
evolved interactions was of order A. That is, low-energy 
phase shifts and binding energies are not distorted if high- 
momentum modes greater than the decoupling scale are 
set to zero (or any arbitrary value) by hand. Therefore, 
it is not surprising that the behavior of e t in terms of A 
or A is qualitatively quite similar. We emphasize that the 
existence of cutoffs where et agrees with the experimen- 
tal value does not imply vanishing 37V forces, as they will 
contribute to other observables. 

The neutron analyzing power A y in n-d elastic scatter- 
ing at neutron lab energy E n = 3 MeV has a maximum 
at the center of mass (cm.) scattering angle 6* cm . = 104 
deg, where the predictions based on realistic interaction 
models underestimate the experimental value by about 
20 %. In Fig. [2] we plot the maximum value of A y as 
a function of RG cutoff A and SRG parameter A. The 
cutoff dependence is quite weak, indicating that this ob- 
servable is not a sensitive probe of short-range force ef- 
fects. The net variation of A y over the range of cutoffs 
is smaller than the discrepancy from experiment of the 
initial interactions, which implies that short-range 37V 
forces are not likely to solve the A y problem. 

The cutoff dependence is even weaker for the n-d 
breakup differential cross section in the space star config- 
uration. We demonstrate that in Fig.[3]for the differential 
cross section close to the center of the space star config- 
uration at E n = 13 MeV; the values measured in two 
different experiments are shown as a reference. These 
flat curves are again an indication that space star cross 
section is not sensitive to short-range physics as already 
found in conventional calculations with different 27V in- 
teractions or by adding a 37V force (ltl [ItJ ■ 
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FIG. 2: (Color online) Neutron analyzing power A y for n-d 
scattering at E n = 3 MeV and # c .m. = 104 deg as function 
of RG cutoff A (left side) and SRG parameter A (right side) . 
The horizontal line at 10 2 A v = 5.86 is the experimental value 
from Ref. fTl. 
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FIG. 3: (Color online) Differential cross section for n-d 
breakup at E n = 13 MeV in the space star configuration 
(50.5°, 50.5°, 120°) at arclength 5" = 6.25 MeV as function of 
RG cutoff A (left side) and SRG parameter A (right side). 
The experimental data are from Ref. (square) and (191 ] 
(circle). 



FIG. 4: (Color online) Total cross section for n- 3 H scattering 
at E„ — 3.5 MeV as function of RG cutoff A (left side) and 
SRG parameter A (right side). The horizontal line at a t = 
2.45 b is the experimental value from Ref. [25| . 
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FIG. 5: (Color online) S- and P-wave contributions to the 
total cross section for n- 3 H scattering at E n = 3.5 MeV. On 
the left side they are shown as functions of RG cutoff A or 
SRG parameter A, while on the right side their correlation 
with the 3 H binding energy is shown. The SRG and RG 
interactions are derived from the AV18 potential. 



III. FOUR-NUCLEON OBSERVABLES 

The results shown in this section are based on the 
solution of the AGS equations 20] in a symmetrized 
form following the technical developments expressed in 
Ref. [U, [53, for all elastic and transfer AN reactions 
below three-body breakup threshold. 

As discussed in Ref. [21(, one of the simplest observ- 
ables in AN scattering is the total n- 3 H cross section at 
that exhibits a resonance around E n ~ 3.5 MeV. This 
peak of the total cross section results from a complicated 
interference between 3 Pj n- 3 H partial waves whose rela- 
tive strength is sensitive to the realistic 2N force one uses. 
While at threshold we find the usual scaling between at 
and e t (a t decreases as |e t | increases), at E n ~ 3.5 MeV 
we observe a breakdown of scaling when we use N3LO 
[2lT ] which is a low-momentum potential when compared 
with the meson-exchange potentials. There N3LO yields 
the largest cross section while not having the lowest \et\- 



Furthermore, in Ref. [24[ it was found that adding the 
Urbana IX 3iV force to AV18 slightly reduces at at the 
peak while more significantly lowering the cross section at 
threshold towards the data as expected through scaling. 

Therefore, in order to investigate the effect of low- 
momentum potentials on a t we plot in Fig. 3] the to- 
tal cross section at the peak versus A (A). In contrast to 
studied 3iV observables, at shows stronger dependence on 
A or A, which is not surprising since the ratio of triples 
to pairs increases. 

In Fig. [5] we split up the total cross section into n- 3 H 
relative S- and P-wave contributions using AV18-based 
Mow k ■ The S'-wave contribution scales well with the 3 H 
binding energy; that scaling is slightly violated for RG 
approach at A < 1.5 fm . In contrast, P- waves show 
no correlation with e t and are responsible for an increase 
of the total cross section at small A (A) values. This is 
consistent with the findings of Ref. [21(. In Fig. [H] we use 
the AV18 potential to show at versus E n for the values of 
A (A) that fit the experimental triton binding energy and 
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FIG. 6: (Color online) Total cross section for n- 3 H scattering 
as function of neutron lab energy for different values of RG 
cutoff A (top) and SRG parameter A (bottom). All results 
are derived from the AV18 potential. The predictions of the 
original AV18 potential (dashed curves) are also shown. The 
experimental data are from Ref. [25( ] . 
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FIG. 7: (Color online) Maximum of the neutron analyzing 
power A v for n- 3 H scattering at E n = 3.5 MeV as function of 
RG cutoff A (left side) and SRG parameter A (right side). 
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for the one that yields deepest binding. While one finds 
that one may describe the total neutron cross section 
over a wide energy range by using A w 1.25 fm -1 in 
RG method or A w 1.8 fin - in the SRG approach that 
also yield reasonable values for et, we emphasize once 
again that these particular values of A (A) do not imply 
vanishing 37V and AN forces, as they will contribute to 
other few- and many-nucleon observables, for example, 
to the ground state energies of light nuclei that do not 
match experiment with those "special" choices of A (A) 

There is a clear correlation between maximum values 
of nucleon analyzing power A y in p- 3 He and n- 3 H scat- 
tering [2l|, [22j|; we therefore study only the latter case. 
Though A y in n-d and n- 3 H scattering are also correlated 
to some extent, their dependence on cutoff is different as 
shown in Fig. [71 it is considerably stronger for n- 3 H. The 
largest increase of A y value at the maximum, by a factor 
1.13 (N3LO) to 1.21 (AV18), is observed around cutoff 
values that yield experimental or deepest binding. How- 
ever, according to Ref. [22], the experimental A y value 
at the maximum for p- 3 He scattering in the same energy 
region is larger than theoretical predictions by a factor 
1.45 (CD Bonn) to 1.55 (AV18). In Fig. [8] we use the 
AV18 potential to show A y versus (9 c m . for the values of 
A (A) that fit the experimental triton binding energy and 
for the one that yields deepest binding. 



FIG. 8: (Color online) Neutron analyzing power A y for n- 3 H 
scattering at E n = 3.5 MeV as function of cm. scattering 
angle for different values of RG cutoff A (left side) and SRG 
parameter A (right side). All results are derived from the 
AV18 potential. The predictions of the original AV18 poten- 
tial (dashed curves) are also shown. 



IV. CONCLUSIONS 

In order to probe the sensitivity of 3N and AN scat- 
tering observables to short-range physics, we used AV18, 
CD Bonn and N3LO based Vi ov/ k potentials that are gen- 
erated through the RG (SRG) method to study their evo- 
lution with the cutoff A (A). Truncating the RG (SRG) 
equations to the two-body level amounts to neglecting 
short-ranged 3N (and higher) forces that are generated 
to preserve exact cutoff independence. Therefore, one 
expects to find residual cutoff dependence in few-body 
observables when only 2N low momentum interactions 
are used. That cutoff dependence may provide a mea- 
sure of the sensitivity of a given observable to omitted 
short-ranged 3N (and higher) forces since the RG evolu- 
tion does not distort the long-ranged forces arising from 
pion exchange, provided A (or A) is well above the pion 
mass, and it is only short-ranged operators that "run" to 
maintain cutoff independence. 

Comparing the results shown in Figs. 12131 with those 
in Figs. 14171 one cannot help noticing that the cutoff dc- 
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pendence of 3iV observables is much weaker than the one 
observed for AN observables. Clearly for the 3 AT observ- 
ables, the cutoff dependence is rather weak, which seems 
to imply that short-ranged 3N forces are not likely to fix 
the two long-standing discrepancies with data mentioned 
above. This is indeed what has been found when the lead- 
ing missing "bare" 3N force, which contains both long- 
and short-ranged operators, is added. Nucleon-deuteron 
A y in elastic scattering and the space star differential 
cross section for breakup barely change by adding a two- 
7r-exchange 3 N force [ill [H, H3| , or an effective 3 N force 
due to the explicit A-isobar excitation [HI |28|, or the 
more recent leading 3 N force from chiral EFT [29| . How- 
ever, there is hope that the subleading long-range 3N 
forces from chiral EFT might be important for the res- 
olution of these problems due to their novel space, spin, 
and isospin structures. 

On the contrary AN scattering observables seem to 
be more sensitive to omitted short-ranged many-body 
forces as demonstrated by the more pronounced depen- 
dence on the cutoff. In n- 3 H scattering at low energy the 
total cross section tr t is dominated by S and P waves 
in the relative n- 3 H motion. The S waves ( 1 S'o, 3 S'i) 
are Pauli repulsive and therefore simply scale with et 
over the whole energy region shown in Fig. [BJ There- 
fore sensitivity to 2N forces comes through the P waves 
( 3 Po, 3 -Pi — 1 Pi, 3 P2) which, in the resonance region, have 



a very complex behavior with the cutoff parameter, lead- 
ing to breaking of scaling with e t . This is consistent with 
the previous findings [2l[ obtained with various 2N po- 
tentials. It also indicates that the at discrepancy may 
be sensitive to missing short or intermediate range 3N 
forces, in contrast to the p- 3 He A y puzzle [22j |. 

From these studies one may conclude that AN scatter- 
ing observables are more sensitive to short range physics 
than the 3N observables where, at low energy, they seem 
to be constrained, to a large extent, by on-shell 2N scat- 
tering and three-particle unitarity, as was expressed l ong 
ago by Brayshaw [301. Recent developments [26|, [27| 
indicate that one needs to fit triton binding energy or 
neutron-deuteron doublet scattering length to constrain 
some other 3N observables that, unlike A y , are sensitive 
to scaling. Nevertheless, this is already fine tuning on top 
of results that are already very close to the experimental 
data. This is not the case for low-energy AN observables. 
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